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Abstract 

We present a study of Higgs liadroproduction through vector-boson fusion 
CN I at the NLO in QCD matched with parton showers. We discuss the matching 

^ ■ systematics affecting this process through a comparison of the aMC@NLO 

• . predictions with the POWHEG and the pure-NLO ones. 
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1. Introduction 

The production of a Standard-Model (SM) Higgs boson (-ff°) through 
the so-called vector-boson-fusion (VBF) mechanism features the second- 
largest cross section among the H^ production channels in hadronic colli- 
sions and, although smaller than the gluon-fusion one by about one order of 
magnitude, it still provides useful complementary informations. After the 
discovery of a SM-Higgs-like particle [l|, |2|, the emphasis is rapidly shift- 
ing towards the determination of its properties, and in this respect VBF 
may play an increasingly important role, owing to its sensitivity to various 
combinations of Higgs couplings [3,], which can be studied by considering 
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different decay channels. However, the very distinctive features of VBF, 
with two jets lying relatively close to the beam line and travelling in op- 
posite directions, render it a challenging case, given that rather severe cuts 
have to be applied in order to reduce backgrounds (among which, in the 
coupling measurement perspective, one may count the contamination due 
to gg — 7- H^). While the typical kinematic regions probed at the LHC 
do not pose problems for perturbative-QCD computations (as shown by the 
behaviour of the rather moderate NLO |j, la, |6|] and NNLO [7|, |8| corrections 
in parton- level results), the presence of two jets in a hadronically-enriched 
environment implies the necessity of using hadron-level simulations such as 
those generated with parton shower Monte Carlo's (PSMC's), in order to 
obtain more realistic predictions. 

It has by now become a rather standard procedure that of matching 
NLO QCD results with PSMC's, by using either the MC@NLO or 



the POWHEG [lO, lUj formalism. Because of the potential importance 
of shower and hadronisation effects and of the good behaviour of NLO 
corrections, VBF appears in fact to be an ideal application for matching 
techniques. However, this has been done so far only in the context of 



the POWHEG approach 12|; in this letter, we amend this by presenting 
MC@NLO results obtained with the fully-automated aMCONLO frame- 
work, and by comparing them extensively with those obtained with the 



code constructed in reference [12| and implemented in the publicly avail- 



able POWHEG-Box framework [13|. The primary motivation for doing 



so is phenomenological. As is known, MC@NLO and POWHEG dif- 
fer by terms of order 0{q^^'^) [ij], i.e. two orders larger than the Horn's; 
furthermore, they differ by logarithmic orders beyond the leading even if 
matched to the same PSMC, owing to the fact that POWHEG generates 
the first emission with own Sudakov form factors, independent of those of 
the PSMQj. While these differences are typically small, consistently with 
their being beyond the nominal accuracy of the calculations, Higgs pro- 
duction in gluon fusion constitutes a striking counter-example, with the 
two approaches yielding significant discrepancies in the Higgs transverse 
momentuno, and in the Higgs-hardest-jet rapidity difference. The latter 
observable in particular, being quite sensitive to the radiation pattern gen- 



""^The latter differences are actually logarithmically leading in the case of an angular- 
ordered PSMC which does not include a vetoed-truncated shower [10|. 
^Before any tuning of the hf act parameter in POWHEG. 



erated by the PSMCo, plus the internal Sudakov in the case of POWHEG, 
could have direct implications for VBF, given the importance of 'extra' ra- 
diation in this process. In general, the differences between the MC@NLO 
and POWHEG results should give one a fair idea of the NLO-matching 
systematics, a topic which, to the best of our knowledge, has not been stud- 
ied in VBF Higgs production. A lesser motivation is technical, and is that 
of validating the aMC@NLO machinery with a further non-trivial process 
on top of those considered so far. 

We remind the reader that aMC@NLO is a generator that implements 
the matching of a generic NLO QCD computation with a PSMC according 
to the MC@NLO formalism; its defining feature is that all ingredients of 
such matching and computation are fully automated. The program is de- 



veloped within the MadGraphS [16| framework and, as such, it does not 



necessitate of any coding by the user, the specification of the process and 
of its basic physics features {e.g. particle masses or phase-space cuts) being 
the only external informations required: the relevant computer codes are 
then generated on-the-fly, and the only practical limitation is represented 
by CPU availability. aMC@NLO is based on different building blocks, 
each devoted to the generation and evaluation of a specific contribution to 



an NLO-matched computation. MadFKS [17| deals with the Born and 



real-emission terms, and in particular it performs, according to the FKS 



prescription [18|, |19[, the subtraction of the infrared singularities that ap- 
pear in the latter matrix elements; moreover, it is also responsible for the 
process-independent generation of the so-called Monte Carlo subtraction 
terms, namely the contributions that prevent any double- counting in the 
MC@NLO cross sections. Finally, MadLoop |20| computes the finite 
part of the virtual contributions, using the OPP |2l| one-loop integrand- 



reduction method and its implementation in CutTools |22 



2. Results 

In this section we present results relevant to the production of a 125 GeV 
Standard- Model Higgs boson through a VBF mechanism at the 8 TeV LHC. 
aMC@NLO includes all interferences between t- and u-channel diagrams, 
such as those occurring for same-flavour quark scattering and for partonic 
channels that can be obtained by the exchange of either a Z° or a W^ boson 
{e.g. ud — )■ H^ud). These interferences, which are kinematically suppressed 



^See reference [15| for a discussion on this point. 
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and de facto negligible, are not included in POWHEG. Furthermore, only 
vertex loop-corrections are considered in both matching schemes, as the 



omitted loops are totally negligible |23| . Electroweak NLO corrections, neg- 



ative and of the order of 5% for this Higgs mass and collider energy 23l. l24| . 
are not included. The Higgs boson is considered as stable. 

Matching with different showers, namely HERWIG6 |25J, HERWIG++ 



26[, and virtuality-ordered Pythia6 [27[ (abbreviated in the following with 
HW6, HWPP, and PY6, respectively), is considered both in aMCONLO 
and in POWHEG, in order to estimate the dependence of physics results 
on the shower model, within the same matching scheme. 



2.1. Setup: parameters and cuts 

Here we list the input settings employed in this computation. The values 
for the Standard-Model parameters follow the prescriptions of the Higgs 
Cross-Section Working Group (HXSWG) Q: 



Mw = 80.398 GeV, 


Tw 


= 2.089 GeV 


Mz = 91.188 GeV, 


Tz- 


= 2.496 GeV 


Gir = 1.166 X 10"^ GeV"^ 







(1) 

Results are obtained by using the MSTW2008NLO PDF set ^, with errors 
estimated at the 68% confidence level. Moreover, renormalisation and fac- 
torisation scales are set equal to the W mass, as suggested by the HXSWG. 
All parton showers are run with their default settings, with the only excep- 
tion of aMC@NL0+PY6, where PARP(67) and PARP(71) are set equal to 
one. Furthermore, no simulation of the underlying event is performed. 

Parton-level events are generated without imposing generation cuts, 
with the exception of a technical cut that requires at least two jets with 
Pt{J) > 2 GeV in the aMC@NLO samples. This cut has been extensively 
checked not to introduce any bias in total rates and differential distributions. 
After shower and hadronisation, typical selection cuts used in experimen- 
tal VBF analyses (called VBF cuts henceforth) are applied: hadrons are 
clustered into jets by using the anti-fc^ algorithm 30| as implemented in 



FastJet [31|], with AR = 0.5. The presence of at least two jets is required, 

with pt{J) > 20 GeV and \y{j)\ < 4.5. Furthermore, the two hardest jets 

{i.e. the two jets with the largest transverse momenta) among those fulfilling 

these criteria are required to have an invariant mass M{ji,J2) > 600 GeV 

and a rapidity separation \Ay{ji,J2)\ > 4. 
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2.2. Differential distributions 

We now present results for various differential distributions. The same 
pattern is adopted for all figures. In the main frame of each figure, the 
three curves that correspond to the aMC@NLO samples are shown: black 
solid for HW6, red dashed for PY6, and blue dot-dashed for HWPP. The 
upper and central insets show, with the same colours and patterns as the 
main frame, the ratios of the aMC@NLO and POWHEG results over 
the fixed-order NLO oneqj, in order to assess the impact of the different 
parton showers and matching schemes on the observables considered. The 
lower insets show the scale (red dashed) and PDF (black solid) uncertainties 
relevant to the aMC@NLO-|-HW6 sample. The scale-variation band is the 
envelope of the results obtained by varying independently the factorisation 
and renormalisation scales in the ranges 

-^ < ^j,, ^p < 2Mw , (2) 



while the PDF errors are computed with the Hessian method [32], as pre- 
scribed by the MSTW set. We remind the reader that the aMCONLO 
Les Houches parton-level event files store additional information sufficient 
to the automatic determination of scale and PDF uncertainties at no extra 



CPU cost, by means of the reweighting technique presented in [33 . 

Part of the differences that will appear in the upper and middle insets is 
due to the different impact of the VBF cuts on the QCD radiation generated 
by the various PSMC's. To better understand this effect, in table [1] we 
quote the ratios of the matched-NLO cross section after VBF cuts to the 
fixed-order NLO one, cr^^^g = 0.388(2) pb. It can be highlighted that these 
ratios are all smaller than one, as typically parton showers tend to spread 
the radiation hardness throughout the phase space, causing slightly more 
events to fail the cuts. On top of this, there is a clear pattern cr^^g > 
^cuTs > ^ott7 both for aMC@NLO and POWHEG. 

In figures [T] to [3] we show the transverse momentum and rapidity of the 
Higgs boson and of the two hardest (tagging) jets. All these observables are 
described with NLO accuracy, as they are non-trivial in their full kinematic 
ranges already at the Born level 0{ag). Therefore, general agreement among 
the two different matching frameworks, as well as among different showers 
is expected. Indeed, all NLO-matched curves are fairly compatible with 



"^Also computed with aMC@NLO and cross-checked against existing resuhs. 

5 





HERWIG6 


PYTHIA6 


HERWIG++ 


aMC@NLO 
POWHEG 


0.93 
0.92 


0.89 
0.86 


0.83 
0.83 



Table 1: Ratios of matched- NLO eross sections to the fixed-order NLO one after VBF 
cuts for aMC@NLO and POWHEG. 

each other once the ratios in table [1] and the theoretical uncertainties in the 
lower insets are taken into account. The comparison with the fixed-order 
NLO prediction, on top of the overall normalisation effect already shown in 
tabledl displays a consistent action of the shower in affecting the jet spectra, 
an effect which is increasingly important as one moves downwards in the 
jet hierarchy {i.e. from the hardest to the softest jets); in fact this trend 
will become even more evident in the case of the third jet (see later). As 
a consequence of the recoil against the shower-enriched jet activity, NLO- 
matched curves display harder and more central Higgs-boson distributions. 
For the observables shown in these figures, PDF and scale uncertainties are 
generally small (typically of the order of ±3% to ±5%), and fairly constant, 
with only mild increases at large transverse momenta. 

Similar conclusions as the ones presented above can be drawn for the 
azimuthal separation between the two tagging jets, displayed in the top 
plot of figure m which also shows excellent shape agreement between fixed- 
order and matched computations. This is reassuring, since this observable 
is particularly sensitive to Higgs-boson quantum numbers as spin and par- 



ity |34j . |35|, |36[, and therefore any theoretical uncertainty is refiected on the 
characterisation of Higgs properties. 

Showering effects are more important for observables such as the in- 
variant mass and the rapidity separation of the two hardest jets. These 
quantities, shown in figure |5l can probe extreme kinematic configurations, 
where the two jets lie at very large rapidities in opposite hemispheres. At 
fixed-order and at the parton level, events with large invariant mass and 
rapidity separation involve partons with energies of up to (9(1 TeV). Such 
energetic partons, when processed through the shower, have a large proba- 
bility to fragment several times, so that the resulting hadron-level jets may 
carry only a small fraction of their energy. This results in the 20%- to 
30%-deficit, with respect to the fixed-order NLO, visible in the matched 
curves at the rightmost edge of these distributions. In the same region, 
the theoretical uncertainties grow up to the level of ±10-15%, especially at 
large rapidity difference. For the PDF's this corresponds to the larger un- 
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certainty at x ~ 1, whereas for the scale uncertainty this can be understood 
as the inadequacy of the choice ^r^f = ^w for such extreme kinemat- 
ics. Indeed we have checked that, by employing a dynamical scale like 
t^R,F = J2iPT,i/'^ (the sum running over final-state partons), the increase in 
the scale-uncertainty band at large rapidity difference is very much reduced, 
to the level of ±5%, with a negligible shift in the central value. 

Observables relevant to the third-hardest jet are more sensitive to the dif- 
ferent matching procedures and to the effects of parton showers, since their 
description at the matrix-element level is only LO. In the lower plot of figure 
mthe exclusive jet-multiplicity is shown. While the 2-jet bin closely follows 
the ratios in table [H the 3-jet bin shows larger differences, with POWHEG 
predicting less events than aMC@NLO. The deficit with respect to the 
fixed-order result in this bin ranges from 15% to 20% for aMC@NLO, 
whereas it is larger, 30% or more, for POWHEG, irrespectively of the 
PSMC employed. This indicates that such an effect is mainly a matching 
systematics rather than being induced by different shower models. It has to 
be stressed, however, that the POWHEG and aMC@NLO predictions are 
still compatible within scale uncertainties, which for the 3-jet bin are about 
±10%, consistently with the LO precision of this observable. From the 4-jet 
bin onwards, the description is completely driven by the leading-logarithmic 
(LL) accuracy of the showers and by the tunes employed, which translates in 
large differences among the various generators. For such jet multiplicities, 
theoretical-uncertainty bands are completely unrepresentative. 

The 3-jet-bin pattern described above determines the normalisation of 
the third-jet transverse-momentum and rapidity distributions, shown in fig- 
ure El These variables can be significantly affected by the different radiation 
produced by the PSMC's. In particular the aMC@NLO results, which are 
quite close to the pure NLO on average, display a ±15% dependence on the 
shower adopted. Conversely, the three POWHEG predictions are slightly 
closer to each other; we reckon that this is a consequence of the hardest emis- 
sion being generated in this formalism by a Sudakov which is independent of 
the actual PSMC employed, as already mentioned in the introduction. The 
POWHEG curves show a 30%-deficit with respect to the pure NLO, and 
are more central and softer than the aMC@NLO ones. Different settings 
in PythiaG have been checked to induce a variation in the results within 
the previously mentioned discrepancy range, which is thus to be considered 
as a genuine measure of the matching systematics affecting these quanti- 
ties. Scale uncertainties are again compatible with the LO nature of these 



observables, of the order of ±10% in the whole rapidity range, and growing 
from ±10% to ±20% with transverse momentum. 

Because of the pecuhar radiation pattern of VBF, which favours QCD 
emissions far from the central-rapidity region, one way to reduce the con- 
tamination due to background processes, as well as from other production 
channels {e.g. gluon fusion), is that of rejecting any event featuring a veto 
jet, namely an extra jet with rapidity lying between those of the two hardest 
jets: 

min{i/(ji), i/(J2)} < Z/(jvcto) < max{y{ji),y{J2)} . (3) 

The predictions for the transverse momentum and rapidity of the veto jet 
are shown in figure [3 The definition in equation ([3]) implies that the more 
central the third jet, the larger the probability that it be the veto jet. Since 
POWHEG predicts a more central third jet with respect to aMC@NLO, 
the veto condition has the effect that the two predictions for the veto jet 
are slightly closer to each other than for the third jet. aMC@NLO yields 
visibly softer and less central distributions, with discrepancies of 20% to 
30% with respect to the pure NLO at small transverse momentum and 
large rapidity. The POWHEG predictions for the transverse momentum 
are softer than aMC@NLO (with the exception of the matching to PY6, 
where shapes are similar), while rapidities are more central, with 20%- to 
30%-discrepancies with respect to the pure NLO at the edges of the spectra. 
As was the case for the third-hardest jet, the observables related to the veto 
jet are described only at LO accuracy, and affected by large uncertainties, 
roughly ±15%. 

3. Conclusions 

We have presented a comparison between aMC@NLO, POWHEG, and 
fixed-order-NLO predictions for VBF Higgs production at the 8 TeV LHC 
This study allows one to asses the NLO-matching systematics affecting this 
process and its various key-observables. Our results can be summarised as 
follows. For observables described with NLO accuracy at the parton level, 
predictions display small theoretical uncertainties, up to ±5%, and show 
a good agreement between the two matching schemes. Differences with 
respect to the pure-NLO predictions result from the action of the shower on 
the jet activity of the events. For quantities described with LO precision, 
theoretical uncertainties are consistently larger, of the order of ±10% to 
±15%, as well as the discrepancy between the two matching prescriptions, 
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with aMC@NLO generally closer to the fixed-order NLO than POWHEG; 
in particular there is a visible effect in the observables related to the third 
jet. Still, all results are largely compatible once theoretical error bands are 
taken into account. 
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Figure 1: Higgs boson transverse- niomcntum (top) and rapidity (bottom) distributions. 
Main frame: aMC@NLO matched with HERWIG6 (black sohd), virtuahty-ordered 
Pythia6 (red dashed) and HERWIG+-I- (blue dot-dashed). Upper (middle) inset: 
ratios of aMCONLO (POWHEG) over the fixed-order NLO, with the same colour 
pattern as the main frame. Lower inset:-! r^cale (red-dashed) and PDF (black solid) 
uncertainties for aMC@NLO-|-HERWIG6.^ee text for further details. 
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Figure 2: Same pattern as in figure [T] for the hardest-jet transverse momentum (top) and 
rapidity (bottom). 
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Figure 3: Same pattern as in figure [T] for the second hardest-jet transverse momentum 
(top) and rapidity (bottom). 
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Figure 4: Same pattern as in figure [T] for tire azimuthal separation of tlie two hardest 
jets (top) and the exclusive jet-multiplicities (bottom). 
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Figure 5: Same pattern as in figure [T] for the rapidity separation (top) and invariant mass 
(bottom) of the two hardest jets. 
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Figure 6: Same pattern as in figure [T] for the third hardest-jet transverse momentum 
(bottom) and rapidity (top). 
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Figure 7: Same pattern as in figure [T] for the veto-jet transverse momentum (top) and 
rapidity (bottom). 
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